The monopole mass in the three-dimensional Georgi— Glashow model 
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We study the three-dimensional Georgi-Glashow model to demonstrate how magnetic monopoles 
can be studied fully non-perturbatively in lattice Monte Carlo simulations, without any assumptions 
about the smoothness of the field configurations. We examine the apparent contradiction between 
the conjectured analytic connection of the 'broken' and 'symmetric' phases, and the interpretation 
of the mass (i.e., the free energy) of the fully quantised 't Hooft-Polyakov monopole as an order 
parameter to distinguish the phases. We use Monte Carlo simulations to measure the monopole free 
energy and its first derivative with respect to the scalar mass. On small volumes we compare this 
to semi-classical predictions for the monopole. On large volumes we show that the free energy is 
screened to zero, signalling the formation of a confining monopole condensate. This screening does 
not allow the monopole mass to be interpreted as an order parameter, resolving the paradox. 
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I. INTRODUCTION 

On the level of classical field equations, the three- 
dimensional Georgi-Glashow model has two phases: 
When the mass parameter of the Higgs field is negative, 
the SU(2) gauge symmetry is broken into U(l), and when 
it is positive the symmetry is unbroken. The phase of the 
system can be determined by a local measurement of, say, 
the scalar field Tr$ 2 , which vanishes in the symmetric 
phase but is non-zero in the broken phase. 

In the broken phase, the field equations have a 
topologically non-trivial solution, the 't Hooft-Polyakov 
monopole Jl], 0, whose energy is concentrated around a 
point-like core. The mass, i.e., the total energy carried 
by a monopole, decreases when the mass parameter ap- 
proaches zero from below, and vanishes in the symmetric 
phase, in the sense that the solution is indistinguishable 
from the trivial vacuum solution. 

In many cases, however, we are more interested in the 
behaviour of the model when fluctuations are taken into 
account. It is immaterial whether the fluctuations are 
thermal fluctuations in a classical field theory or quantum 
fluctuations in a Wick-rotated (2+l)-dimensional quan- 
tum field theory. Both of these systems are described by 
the same partition function, and we shall make no dis- 
tinction between them. Nevertheless, we shall call the 
treatment based on classical field equations "semiclassi- 
cal" even though it is no more accurate in a classical field 
theory at a non-zero temperature than it is in a quantum 
field theory. 

When fluctuations are present, the above simple pic- 
ture changes completely. In particular, the 'symmetric' 
and 'broken' phases are believed to be analytically con- 
nected to each other ||, |], ||, |J. Order parameter can- 
didates that are not gauge invariant, such as ($), vanish 
in both phases, and positive definite observables, such 
as Tr$ 2 mentioned above, are non-zero in both phases. 
It would seem natural that a quantity like the mass of 
a 't Hooft-Polyakov monopole, however, should be pro- 
tected against the effects of the fluctuations by its topol- 



ogy, and that it should therefore serve as an order param- 
eter for the phase transition. If this were the case, the 
phases could not be analytically connected. One example 
of this is the Abelian Higgs model, in which the vortex 
tension indeed acts as an order parameter |?], ||, |j| . 

On the other hand, it is not even obvious that the 
monopole mass can be given a rigorous definition in a 
fluctuating theory, because, in general, one cannot as- 
sume that the field configurations that contribute to the 
partition function are in some sense close to solutions 
of classical field equations. This problem was solved in 
Ref. |pxj| , however, where the monopole mass was defined 
as the increase of the free energy when the total magnetic 
charge of the system is increased by one. Furthermore, it 
was shown how this quantity can be measured in Monte 
Carlo simulations. 

Thus, we have a well defined observable, the monopole 
mass, which could naturally be expected to be zero in the 
symmetric phase and non-zero in the broken phase, and 
still the phases are believed to be analytically connected. 
The purpose of this paper is to explain this apparent 
paradox. 

First, we present a calculation based on a simple di- 
lute monopole gas approximation, which predicts that al- 
though the monopole free energy is indeed non-zero and 
roughly equal to its classical value in a system of inter- 
mediate volume, it decays to zero at exponentially large 
volumes. Therefore, it should actually vanish everywhere 
in the thermodynamic limit. This calculation is very sim- 
ilar to Polyakov's argument Jll], [l^] that the photon has 
an exponentially small mass in the broken phase. 

Second, we measure the monopole free energy directly 
in a Monte Carlo simulation on different volumes using 
the method developed in Ref. |lO| ]. We find that the 
monopole free energy has a volume-independent value in 
a wide range of lattice sizes, which shows that it corre- 
sponds to a localised, point-like object. In agreement 
with the analytical arguments, however, it eventually 
starts to decrease, when the volume is large enough. 

The vanishing of the monopole free energy in the infi- 
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nite volume limit implies that the monopoles condense. 
This leads to confinement of electric charge according 
to the dual superconductor picture and our results 
can therefore be considered as a numerical verification of 
Polyakov's semi-classical argument ]l2| that the Higgs 
phase is confining. In particular, since the monopolc 
free energy vanishes in both phases in the infinite volume 
limit, it does not act as an order parameter, and this re- 
solves the apparent paradox between a smooth crossover 
and the non-analytic behaviour of the monopole mass in 
the semiclassical approximation. 

Within the framework of high-temperature dimen- 
sional reduction |14| ], the three-dimensional Georgi- 
Glashow model is an effective theory for the Yang-Mills 
theory at high temperatures (see, for example, [p_ 5|| and 
references therein). The phase transition of our model, 
however, is not related to the deconfinement phase tran- 
sition of the Yang-Mills theory or QCD. On the other 
hand, our methods can be generalised to four dimensions 
in a straightforward way, and they may therefore be ap- 
plicable also to studying Abelian monopoles in the 
Yang-Mills theory, in particular whether they condense 
at the transition point as has been suggested as a possible 
"mechanism" for confinement p3| . 

Monopole free energies in the Yang-Mills theory have 
been studied before by several groups [fl7t [l8| 19 , |2(J 
using different techniques. In Refs. |l7], |20| fixed bound- 
ary conditions were used to create a monopole, but this 
leads to significant boundary effects. In Refs. [jl8], |T^] 
a monopole creation operator was used, which lets one 
measure not only the mass but also correlation functions 
of the monopole field. With periodic boundary condi- 
tions, however, the operator creates not only a monopole, 
but also an antimonopole somewhere in the system in 
order to satisfy Gauss's law. The advantage of our ap- 
proach is that the system really has a non-zero magnetic 
charge, and because translation invariance is preserved, 
no singularities can arise even near the boundaries of the 
lattice. 

The structure of the paper is as follows. We start by 
discussing the three-dimensional Georgi-Glashow model 
and the lattice definition of its magnetic monopoles in 
Section ||. In Section III, we use semi-classical results to 
motivate our numerical results. We present details of the 
Monte Carlo simulations carried out in Section IV, and 
the results obtained in Section |y|. Finally we discuss our 



findings in Section VI 



The partition function of the theory is formally defined 
as the path integral 



Z 



DQDAi exp 



<TxC 



(2) 



This can be interpreted as a three-dimensional Euclidean 
quantum field theory, or as a classical statistical field 
theory with the Hamiltonian j3H — J d 3 x£. 

The coupling constant, g|, has the dimensions of mass, 
and we can write the parameters of the theory in terms 
of dimcnsionlcss ratios with the coupling constant 



A 



and 



y 



(3) 



(4) 



The notation here reflects the fact that the theory is 
super-renormalisable (in three dimensions), and thus 
only the scalar mass needs a renormalisation countert- 
erm. Even this is only necessary up to the two loop 
level, and its value is known both in the MS scheme |2lf| 
and in lattice regularisation |^2[ ^3|. In Eq. (|4|), m 2 (gf) 
is the MS renormalised mass with renormalisation scale 
M = 3l- 

To study this model in a fully non-perturbative man- 
ner, we formulate the theory in a way that allows nu- 
merical solution by Monte Carlo simulation on a cubic, 
Euclidean lattice consisting of L 3 sites, labelled by a 
triplet of integers x = (x,y,z). The action is given by 
S = ^2g£{x), with the Lagrangian 



C(x) 



II. THE GEORGI-GLASHOW MODEL 

In the continuum, the three-dimensional Georgi- 
Glashow model is defined by the Lagrangian 

L = - Tr (FyFy)+Tr[Dj, $] [A, *]+m 2 Tr $ 2 +A (Tr $ 2 ) 2 , 

(1) 

where $ is in the adjoint representation of the SU(2) 
gauge group, D i = d l + ig z A l and = {ig z )- y \D h D 3 \. 



+ J2{ 2a Tr$ 2 (x) -Tr<5>(x)Ui(x)<P(x + i)Uj (x) 

i 

+m 2 a 3 Tr $ 2 (x) + a 3 A [Tr $ 2 (£)] 2 } , (5) 



where m 2 is the bare lattice mass parameter and (3 — 
4/( a <?i) is the conventional notation for the bare lattice 
gauge coupling. 

We shall treat this lattice theory as an approximation 
to the continuum one, and therefore we parameterise the 
theory in terms of the renormalised continuum couplings 
defined in Eqs. © and (|). We are able to do this, be- 
cause the relationships between the lattice and contin- 
uum couplings are known |22|, but we shall postpone 
discussion of them until SectionjlVj We shall also express 
all quantities in continuum units. 



A. Magnetic monopoles 

It is very well known that in the continuum, the 
field equations have topologically non-trivial solutions, 
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't Hooft-Polyakov monopoles jl], || . They can be charac- 
terised by a non-zero winding number of the Higgs field 
at the spatial infinity, 



N w = 



16ni 



d 2 S k e ljk Tv^> (fijd) (dj&) S Z, (6) 



where $ = $(3> 2 ) -1 / 2 . Although Nw itself is gauge in- 
variant, the integrand is not, and therefore it does not 
have a direct physical interpretation. It can be easily 
seen, however, that Nw actually corresponds to the mag- 
netic charge associated with the residual U(l) gauge in- 
variance. 

To see this, let us define the magnetic field as Q 



n -I 



Tr$F- 



jk 



1 

2i~g 



(7) 



This is a gauge invariant quantity, and agrees with V x A 3 
in the unitary gauge $ = 03. Therefore it is indeed 
the magnetic field associated with the residual U(l) 
symmetry. The corresponding magnetic charge density, 
Pm = V • B, has the following properties: First, because 
Pm is given by a total derivative, the charge inside a given 
volume can be expressed as a surface integral. Therefore 
any local deformation of the fields inside the volume can- 
not change the charge inside the volume. Second, the 
magnetic charge inside a given volume is, in fact, 



47T 

/ M = I d 3 xp M = — N w , 
9 



(8) 



and is therefore quantised in units of 47r/g. These 
two properties imply that the only way the charge in- 
side a volume can be changed is by moving a magnetic 
monopole in or out of the volume. In other words, the 
magnetic charges are topologically stable. 

What is less well known is that these same properties 
are also true for the lattice theory. We can define the 
analogue of Eq. ([?]) as 

Bi = tijkctjk- (9) 
Here ctj k is the lattice U(l) field strength tensor, 

aij = Tru + (x)u l (x)n + (x+i)u j (m)n + (xii{-j)u i l (x+3)n + 

(10) 

and II + = i(l + $). In the continuum limit, ap- 
proaches a 2 Bi. If we define the magnetic charge inside a 
lattice cell as 



Pm{x) 



[Bi (x+i) - Bi {Z) 



(11) 



objects even in a discrete lattice theory, unlike the instan- 
tons of the four-dimensional Yang-Mills theory. Because 
of the quantisation and stability of magnetic charge, it 
makes sense to consider 'microcanonical' partition func- 
tions Zq m which are restricted to configurations with a 
given magnetic charge Qm- The full, 'canonical' parti- 
tion function is then simply 



Z = 



n 



(12) 



Q M =—00 



We define the free energy of a given topological sector by 



In Z, 



(13) 



and the free energy of a monopole as the free energy 
difference of sectors Qm — 1 and 0, 



AF = F x - F 



(14) 



Semiclassically, Zq m = exp[— S(Qm)], where S(Qm) is 
the action of the monopole solution with charge Qm- 
S(Qi) can also be interpreted as the mass of a monopole, 
and with a slight abuse of language we can generalise into 
the fully non-perturbative case by defining the monopole 
'mass' M by 



M = g 2 3 AF 



(15) 



The semiclassical picture would predict that monopoles 
are massive in the broken phase and massless in the sym- 
metric phase. If this were true, the mass would serve as 
an order parameter for the phase transition. 



B. Boundary conditions 

We measure the monopole free energy following the 
method of Ref. |l^] , which for convenience we briefly re- 
view in this Section. 

Our strategy is to work on a finite sized system, and 
impose boundary conditions that force the total magnetic 
charge of the lattice to be either odd or even, whilst pre- 
serving the translation invariance of the system. This is 
important because translation invariance guarantees the 
, of boundary effects. 
Gauss's law rules out periodic boundary conditions, 
as the total charge is constrained to be zero. However, 
translation invariance is preserved by any boundary con- 
ditions that are periodic up to symmetries of the La- 
grangian, and in general they allow a non-zero magnetic 
charge. For instance, 'C-periodic boundary conditions' 



it satisfies the same conditions that guarantee in the con- 
tinuum the topological stability of magnetic monopoles: 
the charge is quantised and can be written as a surface in- 
tegral. These are the same properties that ensure the sta- 
bility of monopoles in the continuum, and thereby mag- 
netic monopoles are well-defined and absolutely stable 



$(n+ij) = -<T 2 $(n)<7 2 = $*(n), U k (n+Lj) = a 2 U k (n)a2 

(16) 

are such that the net magnetic charge can be non-zero, 
but it is constrained to be even 0. We shall refer to 
calculations using such boundary conditions with a sub- 
script '0'. 



Ut 
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Similarly, if the fields are constrained to behave as 

$(n + Lj) = -~<jj&(n)<jj , Uk(n + Lj) = <7jUk(n)<Tj. 

(17) 

on moving around the lattice, the net magnetic charge 
is odd. We term these 'twisted (C-periodic) boundary 
conditions', and denote results so obtained by a subscript 
'1'. It is easy to see that both sets of boundary conditions 
are symmetries of the lattice Lagrangian. 

By a gauge transformation, the twisted boundary 
conditions may be rewritten as (untwisted) C-periodic 
boundary conditions everywhere save at the edges of the 
lattice, where 

U 3 (x,L,L-l) = -U 3 *{x,Q,L-l), 
U x {L-l,L,z) = -U x *(L-l,0,z), 
[/!(!,- l,y,Z) = -Ux*{L-l,y,Q). (18) 

By a suitable redefinition of the fields 



U 3 {x,N,N-l) 
Ur(N-l,N,z) 



-U 3 (x,N,N-l), 
-Ui(N-l,N,z), 
-Ui(N-l,y,N), (19) 



we can express the twisted boundary conditions as a 
theory with C-periodic boundary conditions everywhere, 
but with an additional term in the action that depends 
solely on the gauge fields: 



Zx = J DUiD<Pexp (-S - AS) , 



(20) 



where the change in the action is 



AS = j3 



L-l 



^2 Tr ^23(^,^0,^0) 



x=0 
L-l 



L-l 



^ TrUi 3 (x ,y, zo) + ^2 Tr u M x o, Vo, z$ 1) 
y=o 



z=0 



We emphasise that, because Eq. (|20|) is equivalent to 
Eq. (|l7]) with the translation invariant boundary condi- 
tions, the choice of coordinates [xq, yo, zq) does not affect 
any observable, and in particular, it does not fix the lo- 
cation of the monopole on the lattice. 

In physical terms, AS gives a negative gauge cou- 
pling to three orthogonal stacks of plaquettes which are 
pierced by three mutually intersecting lines on the lat- 
tice. These lines are known in the literature as 't Hooft 
lines JHI. A single, open 't Hooft line creates a pair of 
Dirac monopoles, and has been used to measure their in- 
teraction potential in Refs. @ §7|, |§ |§, f(J. It should 
be noted, however, that Dirac monopoles are rather dif- 
ferent from 't Hooft-Polyakov monopoles. They have 
only half the magnetic charge of the latter, and are sin- 
gular, non-dynamical objects. In our case, the 't Hooft 
lines are closed by the boundary conditions, and therefore 
they do not create any singularities, but a non-singular 
't Hooft-Polyakov monopole. 



The free energy AF (or the 'mass') of a monopole is 
defined by analogy with Eq. ( [l4|) as 



AF = Fi - F n 



hi 



ln(exp(-AS)). (22) 



In the main, however, we shall study the derivative of 
AF with respect to the scalar mass parameter, y, 



1 dAF 

93 dy 



9%V 



(Tr$ 2 )x (Tr* 2 



9l 



9l 



(23) 



where V is the volume of the system. We know that for 
sufficiently large y in the symmetric phase the free energy 
of the monopole will go to zero (at least in the large 
volume limit). If wc see the derivative becoming zero, the 
free energy is at most a constant. In Section IV we also 



measure the free energy at a point in the symmetric phase 
and find it to be consistent with zero. If the derivative 
is zero all over the symmetric phase, it is reasonable to 
assume, then, that the free energy itself is becoming zero. 



III. SEMICLASSICAL EXPECTATIONS 

We now turn our attention to the semi-classical pre- 
dictions for the 't Hooft-Polyakov monopole (see, for in- 
stance, 0), to which we would like to compare our re- 
sults from the fully quantised theory. 

In the broken phase of the theory the scalar field gains 
a vacuum expectation value (VEV) 




(24) 



The semi-classical solution of unit winding number is 
the 't Hooft-Polyakov monopole 0, |), associated with 
an isolated zero of the scalar field. Away from this, the 
scalar field decays towards its VEV, with a characteristic 
length scale 



653 



(25) 



The gauge field simultaneously decays from being SU(2) 
to being asymptotically U(l) with a length scale 



(26) 



« = U 



We thus have a picture where asymptotically the gauge 
fields are Abelian, save within some extended core whose 
size is defined by the above length scales where the gauge 
fields 'unwind' into the full SU(2) gauge manifold. Wc 
shall find that this scenario remains at least qualitatively 
valid when quantum corrections are introduced. 
The mass of this object is, semiclassically, 



T2 = 4tt-/(x) 
9 3 93 



(27) 
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where f(x) is the 't Hooft function. To satisfy the Bo- 
gomolny lower bound on the mass, /(0) = 1. Also, it is 
known numerically (see, for example, |3l]) that for small 
x, f(x) ~ 1 + X. 

The derivative of this mass, as in Eq. (p3|), is 



d M 



2n 



(28) 



Assuming that the monopoles are point-like and non- 
interacting, we can roughly estimate their density to be 



and in general a restriction in the core size by the bound- 
ary will lead to an increase in the (absolute value of the) 
free energy and its derivatives. As a rough estimate, if 
the total non-Abelian flux inside the monopole core is 
roughly constant, then the flux density will vary as the 
inverse of the volume. The total energy of the system 
would then vary as L~ 3 . (The figure 2.5£ is a rough limit 
derived from our results.) 



Intermediate volumes 



9% 



93 



exp 



M 



(29) 



which is suppressed by the exponential of the mass. We 
may also define a mean separation D of the monopoles 
as 



1 



(Dgi) 



9l 



(30) 



When M 3> <?§, there is a hierarchy between D and the 
fundamental length scales £ s and and therefore the 
above assumption of point-like monopoles is valid. 

This is the semiclassical picture for infinite volume. 
What we are interested in is what happens in the quan- 
tised theory of finite volume, and the interplay of the 
system size, aL with the scales above. Particularly, we 
wish to know the fate of the monopole mass on large 
length scales. 



Let us then consider a system that is large enough 
to comfortably accommodate one monopole, but is so 
small that the fluctuations are not likely to create iso- 
lated monopoles (or, more accurately, well separated 
monopole-antimonopole pairs). This is the case when 
£ -C aL -C D. That is, the entropy-action balance is 
dominated by the action cost, which limits us to the mini- 
mum number of monopoles (and antimonopoles) required 
to satisfy the boundary conditions. 

We expect the free energy difference, Eq. (p2|), to be 
that between a system of one monopole and an uncharged 
box. Because a monopole is a localised object, the regions 
far from the monopole core are unaware of the twist in 
the boundary conditions. AF only gets a contribution 
from the monopole core and is therefore independent 
of the volume. In this case, the identification of AF 
with the monopole mass makes sense, and a compari- 
son between the measured values and the semiclassical 
formulae above yields information on the radiative cor- 
rections to the semiclassical monopole. 



A. Small volumes 

Let us first briefly discuss what happens when the vol- 
ume of the system is comparable to, or smaller than, 
the length scales discussed before. The core size of a 
monopole is given by the correlation length £ (we assume 
for sake of argument that £ g and £ s are comparable), and 
therefore if aL < £, there is no room for a monopole in 
the system. If the system is forced by twisted bound- 
ary conditions to contain one monopole, its core will fill 
the lattice and the whole system will be in the confining 
phase. On the other hand, the untwisted system is in 
the Higgs phase. The free energy densities of these two 
phases differ by a certain non-zero amount A/, which is 
essentially the latent heat, and as this is the case in the 
whole volume, we have 

AF w L 3 Af. (31) 

Thus, we can conclude that when aL <C £, AF should 
scale as the volume of the system. 

When the volume is increased, £ < aL < 2.5£, the 
fields start to approach the U(l) of the Higgs vacuum far 
from the monopole core of the twisted system. Nonethe- 
less, the core will be affected by the boundary conditions, 



C. Large volumes 

As the volume is increased such that aL ^> D, the 
entropy gain in introducing well separated monopole- 
antimonopole pairs into the vacuum outweighs the action 
cost and the mean density of topological objects is no 
longer expected to be the minimum commensurate with 
the boundary conditions. The free energy required to in- 
troduce an extra monopole into the system is now less 
than the mass of the single monopole, as we demonstrate 
with a simple model. 



1. The dilute monopole gas 

Following Ref . (32) , where a similar effect was discussed 
in the case of vortices in (2+l)-dimensions, we assume 
that the density of monopoles is low enough, so that the 
probability of finding one in any sub-volume of space is 
independent of whether there are monopoles present else- 
where in the system. In other words, the monopoles are 
assumed to be point-like or that overlap of the cores is 
of vanishing measure. As discussed above, this dilute 
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monopole gas approximation is believed to be valid deep 
in the broken phase. 

The probability of finding n monopoles or anti- 
monopoles (we do not distinguish) in a volume, V, follows 
Poissonian statistics 



p(n;V) = ^ Vr. 



(32) 



We apply this to the volume of the whole lattice V = 
(aL) 3 — g^ 6 (4L/0) 3 . We find different normalisation 
factors for twisted (n G odd, and Afi = sinh(z/oV)) and 
untwisted (n £ even, and A/o = cosh(^oV)) boundary 
conditions. 

The free energy of the system (or its derivative) is ex- 
tensive and the sum of the free energy of the components 
for a dilute gas, and considering the entire system we 
obtain: 



d AF 1 ^ f Po \ 



nGodd 



93/ 



n£even wo 



,.p \ v e~^ v 



(33) 



where po was defined in Eq. @. Eq. @ gives the 
desired plateau for intermediate V, but then decays to 
zero as V — > oo, beginning once V > V c , such that z/ V c = 
1. Note that since ;/oV is simply the typical number of 
monopoles and antimonopoles created by fluctuations, 
this result shows that the monopole free energy decays 
as soon as the fluctuations can create isolated monopoles. 

Crucial in the above calculation is the assumption 
that the monopoles are non-interacting. Although the 
monopoles at least semiclassically have a long range 
Coulomb interaction, we believe this approximation is 
justified, because the interaction is non-confining. Never- 
theless, it is only an approximation, and therefore it must 
be tested in numerical simulations, as we do in Sect. [v| 

We can also see that the above argument would break 
down if we tried to apply it to the four-dimensional case, 
where the monopoles are world lines rather than point- 
like objects. In a Euclidean theory, the action of the 
monopole world line would be proportional to its length, 
and therefore vq would vanish exponentially when the 
limit of infinite time dimension is taken. The same hap- 
pens for vortices in the three-dimensional U(l) theory ||. 
In future work we aim to verify that this is also true for 
vortices in a non-Abelian theory. On the other hand, if 
one of the three dimensions is compact as in the (2+1)- 
dimensional case at a non-zero temperature, vq is finite, 
and again the vortex free energy vanishes in the infinite 
volume limit l|32|. 



2. Confinement 

The prediction of the dilute monopole gas approxi- 
mation that the monopole free energy vanishes in both 
phases in the infinite volume limit is compatible with the 
properties the phase diagram of the theory is believed to 
have. Vanishing free energy means that the monopoles 
condense, and according to the dual superconductor pic- 
ture Jl3| , this gives rise to confinement. 

Indeed, it is known semiclassically that the non-zero 
monopole density gives the photon a non-zero mass even 
in the Higgs phase fllf] , and this leads to confinement. 
Thus it is natural to assume that the Higgs phase is an- 
alytically connected to the confining phase j|, |[. Again, 
this can only be true if the monopole free energy vanishes 
in the Higgs phase, because otherwise it would act as an 
order parameter signalling a transition from the Higgs to 
the confining phase. 

Previous studies || have supported the idea of a 
smooth crossover between the phases, but as they only 
concentrated on local quantities, they cannot be regarded 
as proofs. For instance, in the three-dimensional Abelian 
Higgs model, the phase transition can only be seen in 
practice by measuring non-local observables such as the 
vortex tension or the photon mass j9|. In the present 
case, the predicted non-zero photon mass has not been 
observed in simulations |5|, [| . It is clear from the results 
presented here that the reason for this lies in the very 
large volumes required. 



IV. LATTICE MONTE CARLO SIMULATIONS 

We simulate the Georgi-Glashow model on the lattice 
via Monte Carlo importance sampling of the partition 
functions for both the C-periodic and twisted C-periodic 
boundary conditions. Updates to the lattice were per- 
formed as compound sweeps consisting of one heatbath 
update to the gauge and scalar fields, followed by two 
over-relaxation steps to each. Measurements were made 
once per compound sweep. 

Statistical errors were estimated by jack-knife analysis, 
dividing the data sets into ten bins. For most lattices, 
the bin size was much longer than the autocorrelation 
time of the observables, making them independent. This 
could be seen in an approximate decrease in the statisti- 
cal errors as 1/yN as the number of measurements, N, 
was increased. The only lattice on which this was not 
readily apparent was the p = 4.5, L — 46, where the 
errors did not show such a reduction. This may indi- 
cate that, despite considerable computational effort, the 
ensemble size is still such that the autocorrelation time 
was comparable to the bin size. Error estimates for this 
ensemble should thus be treated as lower bounds. 

To illustrate that the twisted boundary conditions 
( |l7| ) indeed generate a monopole, we show in Fig. |l| 
the isosurfaces of Tr $ 2 and the gauge action density 
53i<_y(l — Trt/y /2) in a typical field configuration at 
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FIG. 1: Isosurfaces Tr$ 2 = 3.89 (green) and £) i<y (l ~ 
TrUij /2) = 0.18 (blue) in a typical field configuration at x — 
0.05, y = 0.45 and /3 = 18. In order to reduce noise, the 
configuration was averaged over 50 subsequent Monte Carlo 
sweeps. 



x = 0.05, y = 0.45 and (3 — 18. The gauge action peaks 
and Tr<I> 2 dips around the same point, exactly as is ex- 
pected to happen near the monopole core. Because of 
thermal fluctuations, the isosurfaces are not spherical. 



A. Observables 

We measure the free energy and its derivative with re- 
spect to the scalar mass. The former is done via Eq. (22). 
In practice this does not work; the importance sampling 
of the theory with untwisted boundary conditions has 
very small overlap with that of the twisted partition func- 
tion. This leads to strong sign fluctuations in AS which 
leads to a poor convergence of its average through Monte 
Carlo simulation. 

Instead, as in Refs. ]|, ||(], p7| , we can introduce a set 
of ensembles defined by a real parameter, s S [0, 1]: 



Z e = J DU t D<P exp (-S - eAS) . 



(34) 



where e = is the untwisted case, and e = 1 represents 
twisted boundary conditions. We then write 



(35) 



where the subscript e indicates that the expectation value 
must be measured using Eq. ([34]). This gives us the ab- 
solute value of AF, but with the cost that we have to 
measure expectation values at non-physical values of e. 



X 


y 


/3 




L 


0.35 


-10, -3, -1, -0.5, 1, 10 


18.0 


0.222 


16 



TABLE I: Lattices used to study the monopole free energy by 
the method of progressive twisting. 



X 


V 


& 




L 


0.35 


-0.124 


4.5 


0.889 


4, 6, 8, 10, 12, 14, 16, 18, 20, 
22, 24, 28, 32, 36, 40, 46 






6.0 


0.667 


4, 6, 8, 10 






9.0 


0.444 


4, 6, 8, 10, 12, 14, 16 






12.0 


0.333 


4, 6, 8, 10, 12, 14, 16, 20 






18.0 


0.222 


4, 6, 8, 10, 12, 14, 16, 20 



TABLE II: Lattices used to study the system size dependence 
of the derivative of the monopole free energy. 



We call this the 'method of progressive twisting'. Calcu- 
lations of the free energy by progressive twisting typically 
used 10,000 to 20,000 measurements for each of 37 values 
of the twisting parameter, e, which are then numerically 
integrated. (For an alternative approach, see Ref. pS|.) 

Alternatively, the derivative of the free energy, 
Eq. (23), may be measured directly, which avoids the 
reweighting problem. We are, however, calculating an 
intensive quantity as the difference of two approximately 
extensive numbers. Maintaining a constant error on 
the former demands increasing accuracy in the latter 
for increasing volume. Even allowing for self-averaging 
and the good scaling properties of the simulation algo- 
rithm, maintaining comparable precision in the free en- 
ergy derivative requires CPU time rising as L 6 . This 
limits the results of this study to L < 46. Calculations 
of the derivative with respect to the scalar mass used be- 
tween 200,000 and 500,000 measurements for each of the 
boundary condition choices. 



B. Lattice parameters 

The physical and lattice parameter values used are 
listed for reference in Tables | and ||. 

In this section we discuss simulations of the SU(2) 
Georgi-Glashow model in three Euclidean dimensions. 
The action for the theory has been given in Eq. (g). In 
addition to the parameters that define our theory in the 
continuum limit, x and y, there are two additional com- 
plications in the lattice theory, being the lattice spacing, 
agf, and the volume, (ag^L) 3 , of the cubic lattice on 
which we perform the simulations. 

Detailed investigations of finite volume effects and scal- 
ing of correlation lengths have been performed for the 
d = 2+1 pure gauge SU (2) and Georgi-Glashow field the- 
ories in jl5|, |3^|. Here we summarise the findings briefly 
for the benefit of non-specialist readers. 

The lattice calculations yield dimensionless results, 
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which may be interpreted as being the physical result 
multiplied by the lattice spacing raised to their naive di- 
mensions, and which we denote via a circumflex accent. 
We remove the dependence on the unknown lattice spac- 
ing by multiplying the result with the appropriate power 
of (3 — 4/ (agf), and therefore it is natural to express the 
results in terms of powers of 53, which has the dimensions 
of (mass) 1//2 . For sufficiently fine lattices, the agreement 
with the continuum limit will be within the statistical er- 
rors of the lattice data, but on coarser lattices there may 
in principle be deviations. The results in Refs. (l^, |3^] 
are indicative of the continuum limit for (3 > 4.5, which 
includes relatively coarse lattices at the lower end of this 
range (as we discuss later). 

The lattice theory in Eq. (||) is parameterised by three 
couplings (m 2 , A, 0). In order to vary the lattice spacing, 
we wish to change (3 whilst maintaining the same contin- 
uum theory [i.e. (x,y)]. This is commonly referred to as 
moving along 'lines of constant physics ' . These trajec- 
tories have been calculated |2^, ^3) in the limit (3 — > 00, 
and they are believed to be valid for lattices finer than 
/?~ 4.5 -5.0: 



A 



y - (4 + 5x) 



3.1759 

47ra 5 2 



16tt 2 



(20a; - 10x 2 ) ( In — + 0.09 

a 93 



+ 11.6x + 8.7 



(36) 



Again, we address the range of applicability in a later 
section. We are primarily interested in testing the idea 
that the 't Hooft-Polyakov monopoles condense. The 
measurement of this is a fine balance. Whilst monopoles 
are topologically stable even if their core is smaller 
than the lattice spacing, it should be much larger than 
that to ensure they resemble the semiclassical 't Hooft 
Polyakov solution. Experience indicates that the corre- 
lation lengths of the gauge and scalar field should be at 
least 2 or 3 lattice spacings. Simultaneously, in order to 
see the screening of the free energy that signals the for- 
mation of the plasma, we require lattices that are (much) 
larger than the mean separation of the monopoles, such 
that it is possible for screening of magnetic charge to 
occur. Given that these two scales may be widely sepa- 
rated, it is not at all clear that we will be able to achieve 
the balance using a lattice size, L, which can be realisti- 
cally simulated on the resources available. 

We can use the known, semi-classical description of 
the monopoles 11 to estimate the parameters needed for 
the lattice. Such estimates are, of course, only expected 
to be accurate up to numerical factors which may be 
important here. Nonetheless, we may hope the results 



are indicative at least, and the exercise gives some insight 
into the possible screening mechanism. 

The monopole density ( p9f ) has a maximum value of 
just under 0.000345. Screening will become apparent 



when the physical volume, g%V 
vqV c = 1; this yields 



> 



4£ c 



L r > 



^0 
3.56/3. 



(4L//3) , is such that 

(37) 
(38) 



If a conservative value of L c = 16 is chosen to allow for 
possible suppression of the monopole density, this indi- 
cates that the gauge coupling is restricted to be (3 < 4.5. 
Our primary interest is in observing the monopole screen- 
ing, so it is not strictly necessary that the perturba- 
tive lines of constant physics still hold on our lattices. 
We would like to maintain some contact with continuum 
physics, however, and thus go no lower than (3 — 4.5. 

Using Eq. (|29|), the maximum monopole density is 
reached for vj gz = 0.421. We are most interested in the 
fate of the monopole mass in the region of the phase dia- 
gram where there is a crossover between the two phases. 
For this reason we select x = 0.35, and thus y — —0.124. 
At this parameter set, the gauge correlation length is, in 
units of the lattice spacing, 



93 



= 2.67, 



and the scalar field correlation length is 



& = T 



= 3.20, 



(39) 



(40) 



which arc both suitably larger than the lattice grid size. 
Finer lattices were used to resolve better the small vol- 
ume behaviour. 



V. RESULTS 



A. Small volumes 



As discussed in Sect. |III A| the free energy AF is ex- 
pected to be proportional to the volume of the system 
when L < £. We studied this in our simulations by 
measuring its y-derivative with couplings x = 0.35 and 
y = —0.124. Obviously, this should behave in the same 
way as the free energy difference itself. The results from 
lattices of different sizes and different lattice spacings are 
plotted in Fig. || as functions of the physical lattice size 
ag\L. At small L the data show very little scaling vio- 
lation. This suggests that we are not seeing a physically 
interesting effect here and supports the idea that the be- 
haviour with L has a simple origin. We show a fit of the 
form —doL dl , where d{o.i} are free parameters. Whilst 
the power law fits well by eye, the precise nature of the 
data makes the fits all quite poor (% 2 /dof > 5). The fit 
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-10 



-30 



-- 



if 



ag 3 =0.222 
ag 3 z =0.333 
ag 3 2 =0.444 
ag 3 2 =0.667 
• ag 3 z =0.889 
power law fit 



x=0.35, y=-0.124 



Sf- 



3 4 



-30 



-40 



P=18, x=0.35 



- c„/sqrt(y -y) 

- y= 

o L=16 
o L=20 
A L=24 



-10 



-6 -4 -2 

y 



FIG. 2: The derivative of the monopole free energy as a func- 
tion of (small) lattice size at x — 0.35, y = —0.124. 



FIG. 4: The derivative of the monopole free energy for inter- 
mediate system sizes for fixed x. Also shown is the a semi- 
classically inspired fit to the data. 



12.0 9.0 



6.0 



4.5 



0.0 0.2 0.4 0.6 0.8 1.0 

(ag 3 ¥ 



FIG. 3: Continuum limit of the derivative of the free energy 
at x — 0.35, y = —0.124 for fixed physical volume, g%L = 5.3. 



shown is to the = 18.0 data only, and gives d\ = 4.2 (5). 
Whilst not precisely 3, this gives qualitative support to 
our simple picture. 

Beyond g\L ~ 3.5 we see different behaviour. The 
derivative now decreases towards a plateau on intermedi- 
ate scales. Whilst this decay may be a power law, we find 
the data insufficient to support a precise fit. The value of 
the plateau does show evidence of a discretisation effect. 
We may attempt to quantify this through a continuum 
extrapolation of the data at g\L = 5.3, admittedly still 
in the transient region, but where we have results for four 
couplings. We show the data in Fig. ||, along with a fit 
assuming only a leading order correction to scaling that is 
quadratic in the lattice spacing. This describes the data 
> 6.0 well (with x 2 / dof = 0.178). Even = 4.5 only 
deviates from this line by 7%, which backs up our pre- 
vious statements on scaling and the applicability of the 
perturbative lines of constant physics (used to maintain 



constant x, y as we varied 0). In addition, this fit sug- 
gests that residual lattice spacing corrections are indeed 
very small at = 18, being around 2% in this case. 

In the region of intermediate volumes, when the 
twisted lattice supports a single monopole, we may at- 
tempt to measure the mass directly, to test the applica- 
bility of the semiclassical results to fully quantised exci- 
tations. We have two methods of approaching this. Less 
prone to statistical uncertainty is to use measurements 
of the derivative of the mass dM/dy over a range in y 
at fixed x. We make a 'mean field' assumption that we 
can describe this data using the formulae of Section III, 
allowing for a shift in the phase transition by the sub- 
stitution y — > y — y c . Typical data, with such a fit, are 
shown in Fig. [|. The mean field assumption fits the data 
well, and from the coefficient cq we may extract a value 
for the radiatively corrected 't Hooft function. We find 
y c to be consistent with zero for x = 0.35. 

Alternatively, we can measure the mass directly by the 
method of progressive twisting for fixed x, y. We show 
such a calculation in Fig. |s|. The dominant error arises 
from the almost complete cancellation of the areas under 
the curve either side of e — 0.5. To illustrate this we plot 
also the same curve rotated through 180°. 

We summarise these estimates of f(x) in Table II], 
and in Fig. ^| where we show a fit to different y as per 
Eq. ( p7| ) . The masses and their derivatives behave much 
as the semiclassical expectations. Similarly the 't Hooft 
function, within the limits of our statistical errors, does 
not appear to differ markedly due to radiative correc- 
tions. There is, however, a considerable variation in the 
data at y = —0.124 as we change 0, and we may worry 
about systematic effects in our results. The first source 
of these is discretisation effects. The majority of our esti- 
mates are for = 18.0, and as we have argued above, the 
residual lattice spacing effects are small here. The varia- 
tion in in the table is also in part due to a correspond- 
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0.2 



0.4 



0.6 



-30 



10 



20 



30 



40 
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FIG. 5: Measuring the monopole free energy by progressive 
twisting of an intermediately-sized system. We show a 180 
degree rotation of the data to highlight the asymmetry. 



X 


method 


y 


M/gl 


/(*) 


/3 


L 


0.05 


deriv 








1.066 (11) 


18.0 


16-20 


0.35 


deriv 








1.257 (14) 


18.0 


16 - 24 


0.35 


prog. 


twist 


-10 


50.8 (14.4) 


1.07 (31) 


18.0 


16 


0.35 


prog. 


twist 


-3 


33.3 (5.5) 


1.28 (22) 


18.0 


16 


0.35 


prog. 


twist 


-1 


18.2 (4.7) 


1.21 (32) 


18.0 


16 


0.35 


prog. 


twist 


-0.5 


13.4 (1.9) 


1.26 (18) 


18.0 


16 


0.35 


prog. 


twist 


1 


1.4 (2.1) 




18.0 


16 


0.35 


prog. 


twist 


10 


1.5 (1.2) 




18.0 


16 


0.35 


prog. 


twist 


-0.124 


2.8 (1.6) 


0.53 (31) 


4.5 


16 


0.35 


prog. 


twist 


-0.124 


5.8 (1.8) 


1.10 (35) 


9.0 


16 


0.35 


prog. 


twist 


-0.124 


6.3 (2.4) 


1.19 (46) 


12.0 


16 



TABLE III: Estimates of the monopole mass and the 't Hooft 
function. 



P=18, x=0.35,L=16 



60 



40 



semiclassical fit 




FIG. 6: The monopole free energy as a function of y for 
intermediately sized system. A semiclassically inspired fit is 
show, giving a value of the 't Hooft function. 



FIG. 7: The derivative of the monopole free energy as a func- 
tion of lattice size at /3 = 4.5, a; = 0.35, y = —0.124. Also 
shown are fits assuming no screening and a dilute gas screen- 
ing picture. 



ing change in the physical volume of the system, and we 
may ask whether all our measurements are for 'plateau' 
masses uncontaminated by the transient small volume ef- 
fects. We believe such biases to be small, especially for 
the y < — 1. As we vary y in Fig. ^ there is a great change 
in the correlation lengths £{s,g} for fixed volume. That 
the different effective volumes considered agree suggests 
we are indeed seeing the intermediate plateau unaffected 
by small L transients. We are thus confident that our 
errors on these estimates of f(x) are accurate. The joint 
fit in Fig. | yields /(0.35) = 1.23(12). 



B. Intermediate and large volumes 

The large to intermediate system size data for the 
derivative of the free energy are shown in Fig. ^. For 
intermediate system size it is clear that the data is well 
represented by a constant independent of the lattice size, 
and we use such a hypothesis: 



1 d M 

9a dy g\ 



Co 



(41) 



where we expect the parameter cq to be po/g 2 - We show 
such fits in Table Our method is to begin with a low 
upper limit for the fitting range, and to then increase 
this, including progressively more data in the fit. The 
X 2 per degree of freedom and Q (if our fitted form is 
the correct model, the probability that our data could 
have arisen as random fluctuations around that model) 
remain (very) acceptable up to L ~ 40. It is clear that 
beyond this the fits become unacceptable: the behaviour 
has changed as a consequence of screening. 

We can attempt to describe the screening by fitting 
over a similar range using a fitting ansatz suggested by 
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-^low 


ihigh 


iVdof 


Co 


X 2 /dof 


Q 


6 


32 


10 


12.431 (58) 


0.686 


0.722 


6 


36 


11 


12.426 (58) 


1.106 


0.353 


6 


40 


12 


12.421 (58) 


1.704 


0.066 


6 


46 


13 


12.419 (58) 


1.813 


0.040 


10 


32 


8 


12.705 (82) 


0.772 


0.628 


10 


36 


9 


12.689 (82) 


1.229 


0.272 


10 


40 


10 


12.677 (82) 


1.873 


0.044 


10 


46 


11 


12.671 (82) 


1.976 


0.027 



TABLE IV: Fitting the derivative of the free energy with an 
unscreened ansatz in Eq. (fi"l|). 



^low 


-^high 


A^dof 


CO 


ci(xl0 5 ) 


X 2 /dof 


Q 


6 


32 


9 


12.499 (44) 


0.151 (93) 


0.938 


0.490 


6 


36 


10 


12.500 (44) 


0.187 (40) 


0.867 


0.564 


6 


40 


11 


12.502 (44) 


0.202 (35) 


0.826 


0.614 


10 


32 


7 


12.445 (83) 


1.36 (106) 


0.816 


0.582 


10 


36 


8 


12.451 (83) 


1.82 (47) 


0.753 


0.675 


10 


40 


9 


12.454 (83) 


1.99 (37) 


0.722 


0.689 


10 


46 


10 


12.449 (83) 


1.79 (35) 


0.773 


0.655 



TABLE V: Fitting the derivative of the free energy with a 
dilute gas screening ansatz in Eq. ((Bgh. 



the dilute gas model: 

1 8 M , n N , e -2ci (a S ^L) 3 

-ott-9 = 4coci(agiLY , 2r , 3 (42) 

where cq is as before, and ci ~ t'o/fff- We show such fits 
over similar ranges in Table [v| For intermediate L the 
fits are similar to those obtained using just a constant 
description. As data from larger systems is included, 
however, we see that the ansatz now remains good. A 
comparison of the two fits is plotted in Fig. [j]. 

We may calculate from ci the mean density of 
monopoles, ^o/fff — 1-3(3) x 10~ 5 , which makes their 
mean separation 

g\D = 42.6 (3.3) (43) 

or, in lattice units at /3 = 4.5, D = 47.8(3.7). From 
this it is clear that we have not got the lattice volume 
necessary to see a complete screening of the free energy 
at L 3> £>■ We cannot therefore completely rule out from 
our data the possibility that AF remains finite even in 
the infinite volume limit. 

As was seen for small system sizes, the plateau values 
at least are heavily influenced by discretisation effects at 
f3 = 4.5. To perform a scaling study of the screening 
mechanism is beyond our current means. Nonetheless, 
for a demonstration of the mechanism such effects are 
immaterial and do not affect the qualitative arguments. 

Note also that no attempt has been made to estimate 
here the systematic errors in the monopole density. To 



do so would require a comparison of different screening 
hypotheses and fit functions, something that the data is, 
unfortunately, not accurate enough to address satisfacto- 
rily. 



VI. CONCLUSIONS 

In this paper, we have used a fully non-perturbative 
technique to measure the free energy of a 't Hooft- 
Polyakov monopole in the three-dimensional Georgi- 
Glashow model. This was achieved by simulating systems 
with two different boundary conditions, both of which are 
periodic up to symmetries of the Lagrangian. This pre- 
serves the lattice translation invariance of the system and 
therefore makes sure there are no boundary effects. 

We found that in the Higgs phase, the free energy 
reached a constant value at intermediate volumes, which 
shows that it is associated with a localised object. 
This is the quantum analogue of the 't Hooft-Polyakov 
monopole. We measure its mass by two different meth- 
ods, and find it compatible with semiclassical expecta- 
tions. 'Mean field' application of the classical relations 
appears successful, and we can make estimates of the 
quantum corrected 't Hooft function. Our best esti- 
mates are /(0.05) = 1.066(11) and /(0.35) = 1.257(14) 
from the derivative of the mass with respect to y, and 
/(0.35) = 1.23(12) by the method of progressive twist- 
ing. These estimates are both self-consistent, and in 
agreement with the classical variation f(x) ~ 1 + x for 
small x indicating that radiative corrections are 

small. 

When the volume increased above a certain critical 
value, however, the free energy started to approach zero. 
This is consistent with an analytical calculation within 
the dilute monopole gas approximation, which predicts 
that the free energy vanishes in the infinite volume limit 
at any values of the couplings. 

In the dual superconductor picture, the vanishing 
monopole free energy implies confinement, and therefore 
our results are numerical evidence for Polyakov's pre- 
diction that the Higgs phase of this theory is confining. 
Furthermore, if the monopole free energy vanishes ev- 
erywhere, it cannot be used as an order parameter, and 
therefore our results strongly support the conjecture that 
the confining and Higgs phases are analytically connected 
to one another. 

Neither, of course, can the monopole mass measured 
from the plateau in the free energy for intermediate sys- 
tem sizes act to distinguish the phases. It is non-zero 
in the deep Higgs phase and zero in the deep symmet- 
ric phase. This plateau does not exist, however, every- 
where in the phase diagram, notably near the transition 
line itself. The mean monopole separation there will be 
comparable to the core size and no plateau would be ob- 
served. Thus the 'mass' is ill-defined and cannot serve 
as an order parameter. 

Our findings suggest a straightforward generalisation 
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to other cases. In a Euclidean formulation in any num- 
ber of dimensions, any point-like topological defect that 
has finite action, will always have a non-zero density at 
any non-zero temperature. This means that these ob- 
jects always have a zero free energy. An extended topo- 
logical defect, such as a string or a domain wall, is, how- 
ever, either a closed loop, surface etc., in which case it 
does not contribute to the global properties of the sys- 
tems, or it has an infinite action. In the latter case, the 
fluctuations cannot generate them, and their free energy 
remains non-zero even in the infinite volume limit. Be- 
cause the free energy can be used as an order parameter, 
this suggests that models with extended topological de- 
fects always have a true phase transition rather than a 



smooth crossover. This question will be studied further 
in a future publication. 
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